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Section  1 


INTRODUCTION 


This  is  the  final  report  under  AFSC  Contract 
F19628-78-C-0073 .  The  subject  of  this  report  is  the 
auroral  E-layer.  The  work  reported  deals  with  its 
modeling  from  a  first,  principles  approach  utilizing  the 
electron  transport  formulation  of  Strickland  et  al . 
(1976)  and  a  time  dependent  chemistry  model  developed 
for  this  program.  A  considerable  amount  of  data  will 
be  presented  in  the  form  of  electron  impact  cross- 
sections,  rate  coefficients  and  parameters  related  to 
optical  emissions.  Results  will  be  in  the  form  of 
electron  and  ion  densities  and  column  emission  rates 
for  a  number  of  UV  features.  Some  of  the  precipitating 
electron  fluxes  considered  for  the  calculations  should 
bear  some  resemblence  to  those  responsible  for  producing 
the  auroral  E-layer  or  equivalently,  the  continuous 
aurora . 


The  report  is  organized  as  follows.  We  begin 
with  a  literature  review  on  the  subjects  of  auroral  elec¬ 
tron  transport  and  chemistry.  The  formulation  applied  to 
this  work  is  then  described.  The  atomic  parameters  used 
in  this  formulation  are  next  discussed.  Included  in  the 
discussion  are  several  tables  and  figures  giving  currently 
applied  values.  We  then  make  comparisons  with  measured 
results  from  two  experiments.  Our  primary  motivation  fo” 
this  was  to  test  the  newiy  developed  chemistry  model. 
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Part  of  the  test  involved  development  of  a  simple  steaciy 
state  model.  Results  from  this  model  are  included  in 
some  of  the  comparisons.  We  finally  direct  our  atten¬ 
tions  to  predictions  which  were  made  to  provide  infoi- 
mation  for  experimentalists  planning  a  program  to  observe 
the  continuous  aurora  in  1981. 
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Section  2 


LITERATURE  REVIEW 


A  brief  review  will  follow  of  theoretical 
studies  carried  out  over  the  past  several  years  con¬ 
cerned  with  electron  transport  and  chemistry  in  the 
ionosphere  under  auroral  electron  bombardment.  Some 
words  regarding  auroral  characteristics  will  first  be 
given  to  provide  an  overview  of  the  problem.  When  an 
auroral  experiment  is  conducted,  it  involves  the  mea¬ 
surement  of  one  or  more  of  the  following  basic  param¬ 
eters  : 

•  non-thermal  electron  (and  possibly 
proton)  distribution  function 

•  densities  and  currents  of  ions  and 
electrons 

•  densities  and  mass  motions  of 
neutral  particles 

•  temperatures 

•  E  field 

•  radiation 

rf 

optical 

X-ray 

The  characteristics  of  these  parameters  are  determined 
by  several  processes.  The  basic  internal  ones  are: 
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•  part  ic  i  e-part  icle  interactions  be¬ 
tween  auroral  electrons  (and  pro¬ 
tons)  and  the  ambient  neutral 
part icles 

•  wave-particle  interactions  between 
auroral  electrons  and  the  ambient 

p  1  asina 

•  chemical  reactions 

•  diffusion  of  the  neutral  and  plasma 
part i c  les 

External  processes  may  also  play  a  role  if  they  generate 
sufficiently  strong  winds,  currents  and  fields, 

A  number  of  other  processes  can  be  identified 
with  respect  to  the  above  list  such  as 

•  ionization 

•  dissociation 

•  excitation 

•  coulomb  scattering 

•  plasma  oscillations  (enhanced  or 
newly  created) 

The-  excitation  of  particles  can  produce  optical  emissions, 
some  of  which  will  be  optically  thick  leading  to  multiple 
photon  scattering.  Coulomb  scattering  of  the  more  ener¬ 
getic  auroral  electrons  will  produce  Bremsstrahlung  radia¬ 
tion  extending  into  the  X-ray  region.  Plasma  oscillations 
can  lead  to  radiation  found  at  rf  frequencies. 

We  will  not  concern  ourselves  in  this  report 
with  wave-particle  effects,  diffusion  effects  on  neutral 
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and  plasma  particles,  or  Bremsstrahl ung  photon  produc¬ 
tion.  We  will  discuss  in  some  detail  particle-particle 
interactions  (excluding  Coulomb  collisions  producing 
Bremsstrahlung)  and  chemical  reactions  as  they  have 
been  incorporated  into  a  Boltzmann  transport  equation 
formulation  being  applied  to  the  AFGL  Auroral  E-Layer 
Program.  A  discussion  will  also  be  given  in  the  appen¬ 
dix  on  ion  diffusion  and  estimates  as  to  when  it  becomes 
important  enough  to  invalidate  local  chemistry  models 
such  as  that  applied  in  this  work.  A  survey  of  selected 
works  on  auroral  electron  transport  and  chemistry  will 
now  follow. 

There  is  an  extensive  literature  on  the  sub¬ 
ject  of  electron  transport  in  the  auroral  ionosphere. 

The  various  approaches  may  be  categorized  as  semi- 
empirical,  range  theoretic,  Fokker-Planck ,  Monte  Carlo 
and  transport  theoretic. 

The  earliest  approach  to  describing  the  trans¬ 
port  properties  of  KeV  auroral  electrons  comes  from 
Rees  (1963)  (a  semi-empirical  approach)  who  applied  an 
energy  dissipation  function  based  on  laboratory  data 
by  Griin  (1957).  The  forms  of  the  functions  obtained  by 
Rees  provide  altitude  profiles  of  the  energy  deposition 
and  ionization  rates  for  monoenerget ic  and  energy  dis¬ 
tributed  sources  with  various  pitch  angle  dependences. 

Walt  et  al.  (1967)  were  the  first  to  provide 
altitude,  energy  and  pitch  angle  information  on  auroral 
electron  fluxes  in  the  KeV  range.  To  do  so  they  ob¬ 
tained  a  numerical  solution  of  a  Fokker-Planck  equation 


which  was  originally  used  to  study  the  properties  of 
electrons  trapped  in  the  radiation  belts  (MacDonald 
and  Walt,  1961).  The  method  assumes  continuous  energy 
loss  and  small  angle  scattering.  Banks  et  al  .  (1974) 
joined  together  Walt's  Fokker-Pl anck  method  and  a  low 
energy  approximate  two-stream  transport  method  by 
Banks  and  Nagy  (1970),  previously  applied  to  photo- 
electron  transport.  The  resulting  equation  was  solved 
numerically . 

Berger  et  al.  (1970,  1974)  chose  to  examine 
auroral  electron  transport  by  applying  Monte  Carlo 
techniques.  Backseatter  yields,  backscatter  spectra 
and  altitude  profiles  of  the  energy  deposition  rate 
are  among  the  transport  quantities  that  were  calculated 
by  this  method.  In  the  first  of  the  two  papers  noted 
above,  information  was  also  given  on  the  lateral 
spreading  of  KeV  electron  beams. 

There  are  at  least  three  transport  models 
utilizing  the  linear  transport  equation  which  have  re¬ 
cently  been  applied  to  auroral  studies.  These  are  the 
models  of  Strickland  et  al .  (1976),  Mantas  (1975)  and 
Stamnes  (1978).  All  three  models  give  a  detailed  de¬ 
scription  oL  elastic  scattering  and  allow  for  discrete 
energy  loss.  Differences  arise  in  the  representation 
of  the  flux  within  the  collision  integral  and  in  the 
method  of  integration  over  depth.  Strickland  et  al . 
(1976;  allow  the  flux  to  vary  quadrat i cal ly  in  S.nE 
and  linearly  in  p  within  any  given  E,p  cell.  The 
quadratic  dependence  was  introduced  because  of  energy 
conservation  problems  for  a  linear  dependence  when 
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treating  energetic  fluxes  above  several  KeV .  The 
integration  over  depth  was  carried  out  by  either  a 
finite  difference  method  (second  order  predictor- 
corrector)  or  an  eigenvalue  method.  The  latter  ap¬ 
proach  was  found  to  be  much  faster  and  more  accurate. 
Mantas  (1975)  allows  the  flux  to  vary  linearly  in 
both  E  and  p  within  any  E,p  cell.  He  also  uses  a 
linear  dependence  within  a  given  z  cell  which  leads 
to  the  standard  finite-difference  expression  for  the 
first  order  z  derivatives.  Stamnes  (1978)  treats  the  u 
dependence  of  the  problem  by  the  discrete  ordinate 
method  and  considers  the  flux  to  be  constant  within 
a  given  E  cell  (commonly  called  the  multigroup  approx¬ 
imation).  For  a  good  discussion  of  the  discrete  ordi¬ 
nate  and  multigroup  methods  see  Davison  (1957).  Like 
Strickland  et  al.  (197G),  Stamnes  (1978)  uses  an  eigen¬ 
value  technique  to  carry  out  the  integration  over 
depth.  All  of  these  transport  methods  lead  to  a  trun¬ 
cated  matrix  equation  which  is  solved  numerically. 

V/e  will  now  briefly  discuss  various  auroral 
ion  chemistry  models.  Most  of  these  do  not  rely  on 
the  transport  results  discussed  above  to  specify  pro¬ 
duction  rates.  Instead,  the  rates  are  estimated  from 
either  optical  (usually  Nri  3914A)  or  ion  composition 
data  with  the  aid  of  relative  strengths  of  electron 
impact  cross  sections.  One  exception  is  the  model 
applied  to  this  program  which  will  be  presented  in  the 
next  section.  Production  rates  are  specified  by  inte¬ 
grating  the  product  of  target  density,  calculated  elec¬ 
tron  flux  and  impact  cross  section  over  electron  energy. 
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It  is  now  well  known  that  the  species  NO  plays 
a  dominant  role  in  auroral  ion  chemistry.  This  knowledge 
came  about  through  the  rocket  experiments  of  Swider  and 
Narcisi  (1970)  and  Donahue  et  al .  (1970).  From  these 
experiments  came  unexpectedly  high  and  low  density 
values  respectively  of  N0+  and  O2  which  led  to  the 
conclusion  that  the  auroral  ionosphere  was  richer  in 
NO  than  previously  thought.  This  was  followed  by  a 
direct  measurement  of  the  auroral  NO  density  by  Zipf 
et  al.  (1970)  which  was  even  higher  than  inferred  from 
the  earlier  ion  measurements.  The  measurement  proved 
to  be  controversial  since  chemistry  modeling  using  con¬ 
ventional  energy  sources  has  not  been  able  to  approach 
the  observed  magnitude  (see,  e.g.  Hyman  et  al.  (1976)). 
Perhaps  the  weakest  link  today  in  the  prediction  of  ion 
densities  in  the  auroral  E-layer  is  the  specification 
of  the  NO  density  unless  given  by  the  experiment  pro¬ 
viding  the  comparison.  A  major  problem  is  its  dependence 
on  energy  deposition  prior  to  the  time  of  interest  which 
follows  from  its  long  lifetime. 

We  will  now  note  some  of  the  modeling  efforts 
over  the  past  decade  and  indicate  how  they  have  dealt 
with  NO.  Swider  and  Narcisi  (1977)  applied  a  steady 
state  model  to  the  analysis  of  ion  data  from  eight 
rocket  experiments.  Since  the  ion  densities  were  given, 
they  were  able  to  infer  the  NO  density.  Their  basic 
conclusion  was  that  modest  enhancements  over  typical 
midLatitude  values  occur  but  far  less  than  reported  by 
Zipf  et  al .  (1970).  Hyman  et  al.  (1976)  carried  out  a 
study  using  a  time  dependent  model  containing  many  ion 
and  neutral  species  with  NO  as  one  of  them.  Production 


rates  were  specified  using  the  electron  transport  model 
of  Strickland  et  al.  (1976).  They  were  also  examining 
the  problem  of  NO  density  enhancements  and  concluded 
that  Zipf's  value  could  not  be  achieved  using  conven¬ 
tional  sources.  Jones  and  Rees  (1973)  reported  results 
from  a  time  dependent  model  which  also  contained  diffu¬ 
sion  of  selected  minor  species.  Production  rates  came 
from  the  deposition  scheme  of  Rees  (1969,  1970,  1975). 
They  presented  ion  composition  and  optical  features  as 
functions  of  time  for  various  model  auroras  noting  that 
auroral  activity  prior  to  the  time  of  interest  can  be 
important  for  some  features  such  as  the  NO  concentra¬ 
tion  . 


One  of  the  most  comprehensive  models  avail¬ 
able  is  that  of  Roble  and  Rees  (1977)  which  joins  to¬ 
gether  Roble's  mid-latitude  F-region  model  (Roble  1975)) 
with  the  above  by  Jones  and  Rees.  The  model  is  claimed 
to  treat  such  effects  as  drift  by  electric  fields,  wind 
motion  and  thermal  energy  and  ion  flow  to  and  from  the 
magnetosphere.  Those  results  reported  are  ion  densi¬ 
ties,  the  Pederson  conductivity  and  the  6300A  and  5577A 
emission  rates  versus  time  without  the  above  complicat¬ 
ing  effects. 

Some  recent  notable  work  has  been  carried  out 
in  connection  with  a  coordinated  rocket-satellite  exper¬ 
iment  reported  by  Rees  et  al.  (1977)  and  Sharp  et  al . 
(1979).  The  model  of  Jones  and  Rees  was  applied  in  the 
analyses  described  in  these  particular  papers.  Predic¬ 
tions  of  ion  densities,  the  electron  temperature  and 
column  emission  rates  for  features  such  as  N2  3914A, 

01  5577A,  01  6300A  and  N2  2P  3371A  were  made  and  compared 
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with  the  data.  Gerard  and  Rusch  (1979)  also  performed 
an  analysis  of  the  ion  density  data  from  the  rocket  ex¬ 
periment  using  the  time  dependent  model  of  Rusch  et  al. 
(1977)  which  allows  for  diffusion.  Production  rates 
are  specified  from  the  incident  electron  flux  by  a 
technique  similar  to  that  developed  by  Rees  (1963).  In 
the  analysis,  the  density  of  NO  was  allowed  to  build  up 
under  various  conditions  and  interact  with  the  other 
species  thereby  providing  long  time  histories  of  key 
quantities  such  as  the  ratio  of  n(N0+)/n(02) .  Good 
overall  agreement  in  ion  densities  was  obtained  by 
allowing  for  substantial  energy  deposition  over  a 
period  of  a  few  hours. 

Val ance-Jones  (1975)  has  also  reported  on  a 
time  dependent  model.  Numerous  ion  and  neutral  species, 
the  electron  temperature,  and  optical  emission  rates  are 
specified.  Production  rates  come  from  the  deposition 
model  of  Rees  (1963)  and  secondary  electron  distribu¬ 
tions  from  the  continuous  energy  loss  description  of 
Rees  (J969).  Diffusion  and  drift  terms  are  not 
included  in  the  coupled  continuity  equations. 


Section  3 


COMPUTATIONAL  MODEL 


In  this  section,  we  will  describe  a  first- 
principles  approach  to  specifying  auroral  electron 
fluxes,  ion  densities  and  column  emission  rates.  The 
block  diagram  in  Figure  1  illustrates  the  steps  and 
inputs  needed.  We  begin  with  an  incident  electron 
flux  versus  energy  and  pitch  angle.  A  Boltzmann 
equation  is  solved  to  obtain  a  steady  state  flux 
throughout  the  important  energy  deposition  region. 

The  needed  inputs  besides  the  incident  flux  are  a  model 
atmosphere  and  electron  impact  cross  sections.  Produc¬ 
tion  rates  for  ions  and  neutrals  are  then  specified  by- 
integrating  the  product  of  density,  flux,  and  cross 
section  over  pitch  angle  and  energy.  These  provide 
the  source  for  a  time  dependent  chemistry  model  which 
currently  gives  altitude  profiles  for  as  many  as  six¬ 
teen  ion  and  neutral  species.  Emission  rates  are  also 
provided  by  the  model  although  most  can  be  specified 
directly  from  electron  impact  on  N2 ,  O2  and  0.  Details 
of  the  above  described  approach  will  now  follow  by  sub¬ 
section  . 


3.1  ELECTRON  TRANSPORT  DESCRIPTION 

The  method  employed  has  been  documented  in 
the  paper  by  Strickland  et  al .  (1976).  A  solution  is 
obtained  to  the  Boltzmann  equation  for  the  following 
conditions: 
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Figure  1  Block  diagram  illustrating  steps  involved  in 
theoretical  modeling  applied  to  the  program. 


•  steady  state 

•  3-D  in  phase  space  (z,  vx,  vB ; 
equivalently  z,E(eV),  indirection 
cosine) 

•  no  external  fields 

•  no  wave-particle  diffusion 

•  discrete  energy  loss 

•  scattering  permitted  through  any 

angle 

•  secondary  electron  production. 

This  equation  may  be  written  as 

-  -KU.t)  <b  +  X. 

cU  f\A 

with  terms  defined  as  follows: 

<t>  electron  flux  in  el/cm2-s-eV-sr 

K  total  inverse  mean  free  path  (IMFP)  (cm-1) 

Kn^  differential  IMFP  for  the  type  of 

particle-particle  interaction  of  the  nth 
impact  species  (N2,  O2  or  0).  Units  are 
cm"  -eV-^ . 

z  altitude 

E  energy 

y  cosine  of  the  pitch  angle  with  respect  to 

the  geomagnetic  field. 

The  differential  IMFP  is  the  product  of  the 
density  and  cross  section: 


r 

, 0) u.' (i) 


niu)^Vu/ 1 


(2) 


Knxu.t',t,e) 


The  cross  section  a  relates  to  the  probability  of  the 
incident  electron  with  energy  E'  producing  an  electron 
of  energy  E  at  angle  G  with  respect  to  the  incident 
direction.  The  outgoing  electron  can  be  either  the 
primary  or  a  secondary  electron.  The  following  inter¬ 
actions  are  modeled: 

•  elastic  scattering 

•  excitation 

•  ionization 

The  cross  section  representation  for  these  processes  is 
given  by  Strickland  et  al .  (1976). 

Equation  (1)  is  replaced  by  a  differential 
matrix  equation  and  then  solved  by  an  eigenvalue  tech¬ 
nique.  As  noted  in  the  earlier  survey,  the  matrix 
form  follows  from  approximating  <|>  in  the  collision 
integral  by  a  function  quadratic  in  JtnE  and  linear 
in  u .  This  functional  form  is  given  over  each  of 
several  (up  to  400)  E,y  cells.  The  code  solving  the 
matrix  equation  contains  tests  for  energy  conservation. 
We  typically  observe  that  the  column  integrated  energy 
deposition  rate  and  the  sum  of  column  energy  rates  for 
lue  many  excitation  and  ionization  processes  are  respec¬ 
tively  within  10%  of  the  net  incident  power. 
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3.2 


SPECIFICATION  OF  PRODUCTION  RATES 


Production  rates  for  ionization,  excitation 
and  dissociation  are  needed  for  the  chemistry  modeling 
and  specification  of  emission  rates.  Various  tech¬ 
niques  have  been  applied  for  their  specification  as 
previously  discussed  in  the  survey  on  chemistry  models. 
Most  have  involved  scaling  techniques  using  relative 
strengths  of  cross  sections  together  with  either  an 
energy  deposition  rate  or  the  N2  ionization  rate  ob¬ 
tained  from  3914A  data.  In  this  work,  the  rates  are 
obtainable  from  the  calculated  electron  flux.  The 
needed  expression  is 

r  * 


,  i 

j 


rvu) 


j  i  ~ 


(3) 


wJL 


where  is  the  spherical  flux  in  el/cm  -s-eV.  The 
sum  allows  for  production  of  an  atomic  species  from 
both  molecules  and  the  parent  atom. 


3.3  CHEMISTRY  MODEL 


A  local  time  dependent  model  has  been  developed 
for  this  program  similar  to  that  reported  by  Hyman  and 
Julienne  (1975).  It  currently  treats  the  following  six¬ 
teen  ion  and  neutral  species:  N0+,  02  ,  0+,  N+ ,  N2, 
0+(2D),  0+(2P),  0£(a4IIu),  NO,  N,  N(2D),  0(1D),  0(1S), 
N2(A3I+),  02(a1Ag)  and  02(b1Eg).  Those  species  whose 
states  have  not  been  identified  are  ground  state  species. 
Atmospheric  model  parameters  held  fixed  during  the  cal¬ 
culations  are  the  N2,  02  and  0  densities  and  the  neutral, 
electron,  ion  and  N2  vibrational  temperatures.  For  a 
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given  incident  electron  flux,  the  electron  and  ion  tem¬ 
perature  profiles  are  estimated  with  the  aid  of  results 
such  as  those  appearing  in  Roble  and  Rees  (1977).  The 
option  is  available  to  either  hold  the  NO  density  fixed 
or  calculate  it.  Finally,  the  source,  which  makes  its 
appearance  in  the  form  of  the  volume  production  rates, 
is  allowed  to  have  time  dependence. 

The  code  which  provides  the  chemistry  descrip¬ 
tion  will  be  referred  to  simply  as  CHEM.  A  feature  of 
CHEM  which  differs  from  most  other  chemistry  codes  is 
its  ability  to  add  or  remove  species  and  reactions 
without  modifying  the  code.  The  chosen  indexing  scheme 
allows  for  such  changes  through  the  input  data. 

A  steady  state  ion  chemistry  code  has  also 
been  developed  for  this  program.  This  was  done  as 
part  of  our  code  validation  effort  on  code  CHEM.  The 
applied  test  involved  running  CHEM  for  long  chemistry 
times  and  comparing  results  with  the  calculated  steady 
state  values.  More  will  be  said  on  such  testing  later 
in  this  report. 

Code  CHEM  solves  a  set  of  coupled  rate  equa¬ 
tions  given  by: 

-  u  eo  (4) 

at  ‘ 


and  its  steady-state  equivalent 
P-(A)  -  (5) 
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where  and  Li  are  the  volume  production  and  loss  rates 
for  the  it*1  species.  A  division  in  species  has  been  made 
between  those  with  slow  and  fast  relative  chemical  life¬ 
times.  In  the  slow  group  are  N0+ ,  0^,  0+ ,  NO,  N,  N(^d) 
and  02(a^Ag).  This  division  follows  the  model  of  Hyman 
and  Julienne  (1975). 

The  processes  modeled  by  the  P  and  L  terms 

include 


electron  impact  (ionization,  dis¬ 
sociation  and  excitation) 

dissociative  recombination 

radiative  recombination 

charge  exchange 

ion-molecule  rearrangement 

neutral  rearrangement 

radiative  de-excitation 

quenching 

The  form  of  either  type  of  term  excluding  electron  im¬ 
pact  (see  Equation  (3))  and  radiative  de-excitation  is 


where  k  is  a  rate  coefficient  which  may  be  temperature 
dependent  and  n{  and  nm  are  the  reactant  species  densi¬ 
ties.  Values  of  important  coefficients  will  be  given 
in  the  next  section. 
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We  include  in  the  Appendix  a  discussion  of 
transport  effects  on  the  ion  densities.  The  study  was 
undertaken  to  assess  the  accuracy  of  the  local  approxi¬ 
mation  in  our  chemistry  model  for  altitudes  near  the 
upper  E-layer  boundary.  Transport  appears  to  have 
little  effect  on  ion  densities  at  these  altitudes. 
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Section  4 


MODEL  INPUT  INFORMATION 


The  atomic  and  atmospheric  data  needed  to  per¬ 
form  the  electron  transport  and  chemistry  calculations 
may  be  divided  into  four  categories.  These  correspond 
to  specification  of: 

•  the  electron  flux 

•  the  production  rates  serving  as 

sources  for  driving  the  chemistry 

•  the  chemistry  determined  ion  and 

neutral  species  densities 

and 

•  the  optical  emission  rates. 

Four  subsections  follow,  one  for  each  of  these 
categories,  which  contain  discussions  and  presentations 
of  some  of  needed  data. 

4.1  DATA  FOR  THE  TRANSPORT  CALCULATIONS 


The  needed  information  is  a  model  atmosphere, 
electron  impact  cross  sections,  grids  in  altitude, 
pitch  angle  and  energy,  and  finally  an  incident  elec¬ 
tron  flux.  Various  model  atmospheres  have  been  applied 
in  this  work.  A  Jacchia  (1977)  model  with  an  exospheric 
temperature  of  1000°K  has  been  used  for  studies  of  a  gen¬ 
eral  nature.  Such  a  study  has  been  one  to  characterize  the 
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electron  density  and  optical  properties  for  the  continu¬ 
ous  aurora  as  reported  in  Section  6.  The  given  model  is 
shown  in  Table  1.  Another  model  applied  in  this  work 
appears  in  Table  2.  The  densities  and  temperature 
shown  come  from  an  auroral  experiment  reported  by  Rees 
et  al.  (1977)  and  Sharp  et  al .  (1979).  Results  using 
this  model  will  be  presented  in  Section  5.  Signifi¬ 
cantly  less  atomic  oxygen  occurs  in  the  latter  model 
which  suggests  there  may  be  important  density  variations 
(from  an  auroral  modeling  perspective)  with  time  in  the 
auroral  ionosphere.  This  must  be  kept  in  mind  for  one 
of  the  problems  of  interest  to  this  program,  namely  the 
specification  of  electron  density  profiles  from  satel¬ 
lite  optical  data. 

Cross  sections  for  electron  impact  on  h'2 ,  O2 
and  0  leading  to  elastic  scattering,  excitation,  disso¬ 
ciation  and  ionization  constitute  most  of  the  needed 
input  data.  Each  such  cross  section  is  input  numerically 
over  an  energy  range  from  1  to  10  eV .  Differential 
information,  as  needed  for  elastic  scattering  and  ioniza¬ 
tion  is  analytically  specified.  We  choose  not  to  present 
the  cross  sections  here  since  they  are  available  in  the 
papers  of  Strickland  et  al.  (1976)  and  Oran  and  Strickland 
(1978)  . 

4.2  CROSS  SECTIONS  FOR  SPECIFYING  CHEMISTRY  SOURCES 

The  initial  species  production  rates  serving 
as  sources  for  the  chemistry  modeling  come  from  electron 
impact  on  N2 ,  O2  and  0.  The  expression  for  their  speci¬ 
fication  was  given  by  Equation  (3).  Table  3  provides  a 
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Table  1  Model  Atmosphere  from  Jacchia  (1977). 


(km) 

Z 

(°  K ) 

T 

a 

(cm-3) 

N2 

(cm  J) 

°2 

(cm  3) 

0 

250 

950 

5 . 5 (  8) 

2 . 9 (  7) 

1 . 5 (  9) 

200 

885 

3 . 2 (  9) 

2.1 (  0) 

4 . 3 (  9) 

150 

664 

3.1 (10) 

2 . 7 (  9) 

1.8(10) 

125 

408 

2.0(11) 

2.1 (10) 

6.5(10) 

110 

242 

1.6(12) 

2.6(11) 

2.3(11) 

100 

194 

9.4(12) 

2.1(12) 

4.6(11) 

90 

190 

5.6(13) 

1.5(13) 

2.4(11) 

30 

210 

3.2(14) 

1.0(14) 

3.2(10) 

Table  2 


Densities  and  temperatures  as  reported  by 
Rees  et  al.  (1977)  and  Sharp  et  al.  (1979). 


Z(km) 

DENSITY 

TEMPERATURE  (°K) 

N2 

°2 

0 

NO 

T 

n 

T 

e 

250 

5.6  (8) 

2.9  (7) 

6.7  (8) 

1.0  (6) 

950 

1900 

200 

4.4  (9) 

4.4  (8) 

1.3  (9) 

8.0  (6) 

886 

1  600 

1  70 

1.3(10) 

1.6  (9) 

2.7  (9) 

1.8  (7) 

790 

1  200 

1  50 

3.5(10) 

4.9  (9) 

5.5  (9) 

7.5  (7) 

660 

920 

130 

1.3(11) 

2.2(10) 

1.4(10) 

3.3  (8) 

472 

590 

1  20 

3.5(11) 

5.7(10) 

2.7(10) 

8.0  (8) 

350 

450 

1 1  0 

1.7(12) 

2.4(11  ) 

6.1(10) 

1.0  (9) 

242 

339 

100 

9.9(12) 

2.0(12) 

1.5(11) 

1.5  (9) 

194 

261 

95 

2.2(13) 

5.5(12) 

1.4(11) 

1.1  (9) 

192 

242 

90 

5.6(13) 

1.5(13) 

8.5(10) 

6.5  (8) 

190 

229 

85 

1.2(14) 

4.0(13) 

4.0(10) 

2.2  (8) 

198 

2  00 

80 

3.2(14) 

1.0(14) 

1 .0(10) 

1.0  (8) 

210 

210 
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list  of  the  species  produced  at  least  in  part  by  elec¬ 
tron  impact  and  their  parent  species  in  this  process. 

The  assumed  cross  sections  for  each  impact  process  are 
given  in  Figures  2-4. 

The  cross  sections  for  species  without  spec¬ 
troscopic  notation  refer  to  effective  production  in 
their  ground  states.  Effective  ground  state  produc¬ 
tion  for  is  assumed  to  be  given  by  the  cross  sec- 

“  + 

tion  for  total  N2  production. 

The  information  shown  in  Figures  2-4  is  based 
on  data  and  estimates  appearing  in  Rapp  and  E.  Golden 
(1965),  Mantas  (1973),  Myers  and  Schoonover  (1975),  and 
Oran  and  Strickland  (1978). 

4.3  RATE  COEFFICIENTS 

Sixteen  species  are  currently  modeled  as 
identified  in  Section  3.3.  Tables  4  and  5  give  the 
dominant  reactions  and  rate  coefficients  for  those 
bearing  directly  on  the  ion  concentrations.  The  species 
O(lD),  O(lS),  and  (^(b^I)  have  not  been  included.  The 
information  provided  has  come  from  Roble  and  Rees  (1977), 
Hyman  and  Julienne  (1975),  Gerard  and  Rusch  (1979), 

Swider  and  Narcisi  (1977),  and  Wolfsy  and  McElroy  (1977). 
Most  of  the  coefficients  may  be  found  in  the  first  two  of 
the  above  cited  references. 
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Table  4  Reactions  and  Rate  Coefficients  for 
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See  NO  Production 
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Table  5  Reactions  and  Rate  Coefficients  for 
NO,  N,  N(2D),  and  N0(A3E) 
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Table  5  Reactions  and  Rate  Coefficients  for 
NO,  N,  N(2D),  and  N„(A3Z) 


4.4  OPTICAL  EMISSION  DATA 

We  currently  calculate  volume  and  column  emission 
rates  for  the  following  features:  Ng  IK  3914A, 

N0  LBH  1325  A,  LBH  1384  A,  N  2P  3371  A,  NI  1200  A, 

NI  1134  A,  KI  1493  A,  01  1304  A,  01  1356  A,  01  2972  A, 
and  01  5577  A.  Multiple  scattering  effects  are  important 
Cor  NI  1200  A,  NI  1134  A,  and  01  1304  A,  and  for  side 
viewing  situations,  there  is  also  some  effect  on  LBH 
bands  and  01  1356  A.  We  assume  that  only  the  01  5577  A 
feature  is  affected  by  chemistry. 

Tables  6-8  contain  information  on  electron 
impact  and  photon  absorption  processes  important  to  the 
modeling  of  the  above  features.  Figures  5-7  give  the 
corresponding  electron  impact  cross  sections.  Not  all 
of  the  features  noted  in  these  tables  and  figures  will 
be  considered  when  presenting  results  below'.  Excluded 
will  be  those  noted  above  strongly  affected  by  multiple 
scattered  and  01  5577  A. 


TABLE  7.  ATOMIC  NITROGEN  EMISSION  INFORMATION  FOR 
ELECTRON  IMPACT  EXCITATION 


TABLE  8.  ATOMIC  OXYGEN  EMISSION  INFORMATION  FOR 
ELECTRON  IMPACT  EXCITATION 
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9.  Stone  and  Zipf  (1974)  10.  Mumma  and  Zipf  (1971) 

1.  V,-]ls  ut  al.  (1971)  12.  Henry  et  al  .  (1969) 
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Section  5 


COMPARISONS  WITH  EXPERIMENT 

In  this  section,  we  compare  electron  densities 
and  other  selected  quantities  obtained  from  our  time 
dependent  chemistry  model  with  those  from  two  exper iments , 
In  the  first,  ttie  measured  density  comes  from  data 
obtained  with  the  Chatinika  incoherent  radar  backscatter 
facility  (Vondrak,  private  communication).  The  needed 
precipitating  electron  flux  data  came  from  a  rocket 
experiment  by  Evans  and  Moore  (1979).  Both  experiments 
recorded  their  measurements  in  the  same  vicinity  and  at 
the  same  time.  The  auroral  activity  was  compatible  with 
that  characterizing  the  continuous  aurora. 

The  second  measured  electron  density  comes  from 
a  coordinated  rocket-satellite  experiment  reported  by 
Rees  et.  al .  (1977)  and  Sharp  et.  al .  (1979).  Several 
parameters  were  measured  including  ion  and  neutral 
densities,  the  precipitating  electron  flux,  and  intensities 
of  optical  features.  One  would  not  classify  the  observed 
aurora  as  a  continuous  aurora  although  conditions  were 
exceptionally  stable  considering  how  energetic  the 
electron  flux  was  observed  to  be. 

COMPARISON  1 

The  'i ecipitating  electron  flux  used  in  the 
transport  calculations  is  shown  in  Figure  8  and  comes  from 
Evans  and  Moore  (1979).  The  distribution  above  'w  1000  eV 
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1 


is  similar  to  a  Maxwellian  distribution  (characteristic 

energy  of  'v  2000  eV )  which  is  typical  for  the  continuous 

aurora.  The  energy  content  of  the  given  flux  is  %  2  ergs/ 
2 

cm  -  s  for  isotropy  over  the  downward  hemisphere  which  is 

assumed  for  the  calculations.  The*  model  atmosphere  in 

Table  1  has  been  used  in  the  calculations.  Figure  9  gives 

the  comparison  in  electron  density  profiles.  A  family  of 

curves  is  presented  for  the  calculations,  each  curve 

referring  to  a  different  integration  time.  The  maximum 

7  -3 

NO  density  from  the  calculations  is  %  5  x  10  cm  .  This 
is  essentially  the  input  value  since  NO  changed  little 
over  the  time  interval  considered.  We  see  that  the  model 
is  able  to  reasonably  reproduce  the  shape  of  the  data,  but 
predicts  a  maximum  density  value  half  as  large.  This  has 
led  us  to  examine  the  applied  chemistry  model  by  comparing 
its  results  under  similar  conditions  with  those  from  a 
steady  state  model.  A  discussion  of  the  comparison  will 
follow  shortly.  Based  on  this  comparison,  excellent  energy 
conservation  in  the  electron  transport  results,  and  the 
general  insensitivity  of  the  derived  electron  density  to 
the  range  of  uncertainty  in  the  chemistry  input  parameters, 
we  would  suggest  that  most  of  the  above  discrepancy  comes 
from  one  or  more  of  the  following  sources  or  conditions: 

•  measured  electron  flux 

•  radar  data  analysis 

•  different  spatial  resolutions  of  the 

experiments . 

The  last  item  is  probably  a  strong  candidate  for 
the  source  of  discrepancy. 

We  now  turn  to  the  discussion  of  the  steady  state 
model.  The  case  to  be  addressed  in  comparison  of  results 
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DASHED  CURVE  GIVES  THE  RADAR  MEASURE 


between  this  and  the  time  dependent  model  is  the  one  having 
the  incident  electron  flux  shown  in  Figure  8.  We  choose  to 
consider  the  ions  which  play  the  dominant  role  in  auroral 
chemistry,  namely  N0+ ,  0+ ,  0+  and  N*.  The  ion  is  itself 
not  an  important  component  to  the  total  ion  density  but  is 
important  as  a  source  of  ionization  for  driving  the  chem¬ 
istry.  We  consider  the  following  equations  for  the  above 
species : 

N0+ 

K]Hki +  KII°K  -  lN0+l  H“i  <s> 


N2 1  [  °2 1  k4  +  (°1  1  °2  1  k5  +  P0+ 


Kilo  |  *2  *  I  °2 )  |N0  I  "l 


0 


(7) 


[nJ]  [o]  k3  +  P0+  =  [o+] |02]  k5 


N 


+ 

2 


PN+  = 


Nr 


0r 

4 


411° 


(kr 


V 


(9) 


Elect ron  density 


(10) 
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These  enable  immediate  determination  of  N*  and  0+ , 
namely  by  rewriting  equations  (9)  and  (10)  as 


Kl  *  W  |K1  k4  +  I°]  <k2  *  k3>)  (11) 


and 


0 


(Kl  l°l  k3  +  v)/|l° 


2  J  k5 


(12) 


With  these  concentrations  specified  and  using  equation  (10) 
equations  (6)  and  (7)  may  be  joined  to  give  a  cubic  equa¬ 
tion  (a  quartic  which  reduces  to  a  cubic)  in  either  N0+ 
or  0*.  Only  one  root  provides  a  reasonable  solution  which 
may  then  be  used  to  find  the  remaining  concentration  by 
solving  a  quadratic  equation  in  this  variable. 


The  inputs  to  the  model  are  the  densities  of 
N2.  09,  0,  and  NO,  the  production  rates  of  N* ,  O* ,  and  0+ , 
and  the  various  rate  coefficients.  The  chosen  densities 


Q)  and  O  are  the  same  used  in  the  time  dependent 


of  .\:2,  0 

calculations.  Three  sets  of  NO  densities  have  been  con- 
si  tiered.  One  of  these  appears  in  Table  9  along  with  the 
assumed  electron  temperature  profile  and  density  profiles 
of  N2,  02  and  0.  A  second  set  has  densities  that  are  ten 
times  the  magnitude  of  those  in  the  table.  A  third  set 
lias  zero  valued  densities.  The  needed  rate  coefficients 
are  among  the  larger  set  appearing  in  Table  4.  Table  10 
i ves  the  production  rates  which  are  the  same  ones  used 
above  in  the  time  dependent  calculations. 


Ion  densities  have  been  obtained  from  both  the 


time-dependent  and  t. ime- i ndependent  formulations  for  the 
above  input  data.  The  time-dependent  scheme  was  simplified 
to  treat  only  the  four  ion  species  and  consider  only  the 


TABLE 

9.  MODEL 

ATMOSPHERE 

(km) 

(  °  K  ) 

(  °  K  ) 

(cm  3) 

(cm  3) 

(cm-3) 

(cm  3 ) 

Z 

T 

T 

No 

0_ 

0 

NO 

a 

e 

2 

2 

250 

950 

2900 

5 . 5  (  8) 

2 . 9 (  7) 

1  .  5 (  9) 

1.2(6) 

200 

885 

2200 

3 . 2 (  9) 

2.1 (  8) 

4 . 3 (  9) 

7.0(6) 

150 

664 

626 

3.1  (10) 

2 . 7  (  9) 

1.8(10) 

4.2(7) 

125 

408 

320 

2.0(11  ) 

2.1 (10) 

6.5(10) 

7.0(7) 

110 

242 

240 

1 .6(12) 

2.6(11  ) 

2.3(11) 

7.2(7) 

100 

194 

210 

9.4(12) 

2.1(12) 

4.6(11  ) 

6.0(7) 
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TABLE  10.  ION  PRODUCTION  RATES  BASED  ON  THE 
ELECTRON  FLUX  DATA  OF  EVANS  AND 
MOORE  (1979). 


(km) 

1 

P  + 

N2 

-1 

(cm  -S  ) 

S 

V 

250 

3.7(1) 

5 . 2  (  -1 ) 

3.1(1) 

200 

1.4(2) 

2 . 5  (  0) 

5.9(1) 

150 

7.9(2) 

1 . 9  (  1) 

1.6(2) 

125 

2.8(3) 

8.0(  1) 

3.6(2) 

110 

5.6(3) 

2 . 5 (  2) 

5.0(2) 

100 

3.9(3) 

2 . 4 (  2) 

3.0(2) 
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reactions  and  corresponding  rates  treated  by  the  time- 
independent  scheme.  Furthermore,  integration  was  carried 
out  to  100  min  for  the  time-dependent  case  to  insure  that 
steady  state  had  been  achieved  over  the  chosen  altitude 
range  from  250  to  100  km.  Integration  times  to  reach 
steady  state  are  actually  much  shorter  than  this  at  the 
lower  altitudes. 

Figure  10  shows  the  calculated  density  profiles 
for  N0+ ,  Og  and  0+  for  the  three  chosen  sets  of  NO  densit¬ 
ies.  The  solid  curves  give  the  time- i nt egra ted  results 
while  the  symbols  give  the  time-independent  results.  The 
fourth  density,  N+ ,  is  about  two  orders  of  magnitude  below 
the  summed  value  and  for  this  reason,  we  have  chosen  not 
to  include  it  here.  The  densities  above  about  200  km  are 
suspect  due  to  the  absence  of  the  effect  of  ion  diffusion 
in  the  formulations.  We  chose  to  carry  out  the  calcula¬ 
tions  to  250  km  just  to  see  how  the  concentrations  behaved 
in  the  local  approximation  where  0+  becomes  the  dominani 
ion . 


We  observe  that  there  is  excellent  agreement 
between  the  two  formulations  which  demonstrates  that 
the  time-dependent  chemistry  code  is  properly  function¬ 
ing  for  its  given  set  of  inputs.  We  further  see  that 
variations  in  NO  weakly  effect  the  total  ion  or  electron 
density  but  strongly  affect  the  0+  concentration.  This 
is  well  known  and  is  simply  explained  by  the  fact  that 
NO  provides  an  effective  means  of  transferring  the  charge 
from  Og  to  NO,  thus  enhancing  NO  while  depleting  Og . 
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FIGURE  10  EQUILIBRIUM  ION  DENSITIES  BASED  ON  THE  ELECTRON 
FLUX  SPECTRUM  IN  FIGURE  8.  EACH  PANEL  REFERS  TO 
A  DIFFERENT  SET  OF  NO  VALUES.  THE  RESULTS  IN 
THE  FIRST  PANEL  ARE  FOR  NO  =  0.  THE  SIMPLIFIED 
CHEMISTRY  MODEL  DESCRIBED  IN  THE  TEXT  HAS  BEEN 
USED.  THE  SOLID  CURVES  GIVE  THE  TIME  INTEGRATED 
VALUES  AT  100  MIN  WHILE  THE  SYMBOLS  GIVE  THE 
RESULTS  FROM  THE  STEADY  STATE  MODEL. 


COMPARISON  2 


An  extensive  set  of  data  has  been  reported  in  the 
literature  for  a  very  stable  aurora  generated  by  an  ener¬ 
getically  hard  electron  spectrum  (Rees  et  al.  (1977), 

Sharp  et  al .  (1979),  and  Sharp  and  Torr  (1979)).  Although 
the  spectrum  of  the  precipitating  electrons  is  much  harder 
than  those  to  be  found  in  the  continuous  aurora,  we  have 
chosen  to  make  comparisons  with  the  data  lor  ion  densities 
and  optical  emissions  since  the  experiment  appears  to  be 
well  enough  controlled  and  conditions  stable  enough  to 
effectively  test  our  chemistry  model. 

As  reported  by  Rees  et  al .  (1977),  very  little 
temporal  variation  in  the  aurora  occurred  over  a  20  minute 
period  based  on  ground-based  photometric  measurements. 
During  this  time  ('v  2200  LT  on  March  20,  1974),  a  rocket 
was  launched  from  Ft.  Churchill,  Manatoba  into  the  display 
while  the  Atmospheric  Explorer  C  satellite  passed  nearly 
directly  overhead  recording  densities,  auroral  electrons, 
and  optical  emissions.  Similar  measurements  were  made 
from  the  rocket. 

Table  2  gives  neutral  densities  and  temperatures 

reported  in  the  above  papers.  These  have  been  used  in  the 

calculations  discussed  below.  The  incident  flux  for  the 

calculations  is  shown  in  Figure  11  based  on  the  AE-C  data 

appearing  in  Figure  2  of  Rees  et  al .  (1977).  The  energy 

2 

content  of  this  flux  is  ^8.8  ergs/cm  -s.  The  solid  portion 
of  the  curve  represents  a  smoothed  version  of  the  data 
while  the  dashed  portions  are  extrapolations  provided  by  us. 
The  electron  and  ion  densities  as  measured  from  the  rocket 
are  shown  in  Figure  12.  They  come  directly  from  Figure  2 
of  Sharp  et  al .  (1979).  Altitude  profiles  of  optica) 
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ELECTRON  FLUX  (el / CM  -S-eV-SV) 


1 


ENERGY 


FIGURE  11  INCIDENT  ELECTRON  FLUX  USED  TO  PERFORM  THE 
CALCULATIONS.  THE  SOLID  PART  OF  THE  CURVE 
IS  A  SMOOTHED  REPRESENTATION  OF  DATA 
APPEARING  IN  REES  ET  AL  (1977).  THE  DASHED 
PORTIONS  COME  FROM  EXTRAPOLATIONS  BY  US. 
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ALTITUDE  (KM) 


FIGURE  12  MEASURED  ION  DENSITIES  APPEARING  IN  SHARP  ET  AL.  v'1979) 


emissions  have  been  reported  by  Sharp  et  al .  (1979)  for 
the  following  features:  IN  3914A,  VK  3220A,  2P 

3371 A ,  NI  5200A ,  01  5577A ,  and  01  6300A.  Figure  13  shows 
the  measured  profiles  for  N*  3914A  and  2P  3371A.  The 
atomic  features  have  been  excluded  since  we  are  concen¬ 
trating  our  attention  on  the  UV  to  XUV  spectral  region. 

We  further  exclude  the  VK  feature  since  we  are  not  current¬ 
ly  considering  VK  bands  in  our  modeling.  The  solid  curves 
in  Figure  13  are  calculated  results  to  be  discussed  shortly. 

Ion  production  rates  based  on  the  incident  flux 
in  Figure  11  and  used  as  inputs  for  the  chemistry  calcul¬ 
ations  are  shown  in  Figure  14.  The  resulting  important 
ion  densities  are  shown  in  Figure  15  along  with  the  measured 

9  _a 

values.  The  maximum  NO  density  is  1.5  x  10  cm  (see 
Table  2)  and  20  minutes  was  allowed  for  the  integral ion  t  ime 
We  regard  the  agreement  shown  as  satisfactory  considering 
possible  uncertaini tes  in  both  the  experimental  and  model 
parameters.  It  does  appear  that  significantly  less  struc¬ 
ture  occurs  in  the  measured  profile  although  we  would 
expect  a  sharp  rise  had  it  extended  to  lower  altitudes. 

We  are  not  able  to  reproduce  the  double  peaks  in  the  total 
ion  (or  equivalent  electron)  concentration.  This  may  result 
from  variations  recorded  by  the  experiment  as  the  rocket 
moved  across  field  lines.  Time  variations  seem  less  likely 
since  the  aurora  was  particularity  stable  during  the  flight. 

Returning  to  Figure  13,  we  present  our  calcula¬ 
tions  of  zenith  viewing  intensities  for  various  UV  features. 
Good  agreement  with  the  data  is  achieved  for  maximum  inten¬ 
sity.  The  calculated  profiles  fall  off  more  rapidly  with 
increasing  altitude,  which  is  likely  due  to  a  more  rapid 
decrease  in  the  model  N£  density  or  too  few  incident 
electrons  in  the  low  keV  range  compared  to  actual  conditions. 
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CALCULATED  VOLUME  PRODUCTION  RATES  OF  IONS  USED  AS 
THE  CHEMISTRY  CALCULATIONS. 


DENSITY  (CM  > 


The  fall-offs  at  low  altitudes  for  01  1356A  and  the  N,^ 

LBH  bands  are  due  to  photon  absorption  by  0 Absorption 
is  stronger  at  1384A  compared  to  1325  A  which  leads  to 
significant  differences  in  the  two  band  altitude  profiles. 
This  effect  can  be  useful  in  attempting  to  estimate  the 
spectral  characteristics  of  precipitating  electrons  from 
optical  data. 
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Section  6 


PREDICTIONS  FOR  THE  CONTINUOUS  AURORA 


In  this  section,  production  rates,  ion  densi¬ 
ties,  the  electron  density,  and  column  emission  rates 
will  be  presented  for  incident  electron  fluxes  given  by 
Maxwellians  with  characteristic  energies  of  500,  1000, 
2500,  and  5000  eV .  The  power  content  of  the  incident 

p 

electrons  for  all  cases  has  been  chosen  as  1  erg/cm  -s. 
The  given  distributions  have  been  selected  to  simulate 
the  types  characteristic  of  the  continuous  aurora. 

In  performing  a  set  of  calculations  as  follows 
to  provide  predictions  to  those  who  conduct  auroral  mea¬ 
surements,  a  number  of  parameters  need  be  kept  in  mind. 
Among  these  are  the  densities  of  0,  Op,  and  NO  and  the 
electron  temperature.  Such  quantities  can  vary  signif¬ 
icantly  in  the  auroral  ionosphere  and  do  have  measurable 
effects  on  ion  densities  and  optical  emissions.  Their 
variations  are  less  important  in  determining  the  electron 
density.  A  single  set  of  atmospheric  parameters  has  been 
used  for  results  to  follow.  The  model  atmosphere  applied 
is  that  given  in  Table  9.  The  NO  density  was  allowed  to 
vary  but  did  not  deviate  significantly  from  the  chosen 
initial  values  for  the  integration  times  allowed  and 
modest  energy  content  of  the  source  electrons.  One  can 
estimate  the  effect  of  varying  the  NO  density  by  examin¬ 
ing  the  results  in  Figure  10.  The  given  model  atmosphere 
is  probably  rich  in  0  based  on  measured  auroral  values 
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such  as  those  appearing  in  Table  2.  The  strongest  effect 
from  this  will  be  in  enhancement  of  01  emissions  such  as 
1356  A. 


The  incident  electron  fluxes  for  the  four  cases 
to  be  considered  are  shown  in  Figure  16.  Isotropy  is  as¬ 
sumed  over  the  downward  hemisphere  at  the  upper  boundary. 
Figure  17  gives  the  total  ion  production  rates  corres¬ 
ponding  to  the  four  Maxwellian  electron  distributions. 

An  indication  of  individual  ion  species  contributions 
may  be  obtained  from  Figure  14. 

The  calculated  electron  densities  appear  in 
Figure  18  for  a  20  minute  integration  time  which  effec¬ 
tively  gives  steady  state  values  except  near  the  upper 
boundary.  A  significant  change  in  shape  occurs  in 
going  from  the  2500  eV  to  the  1000  eV  Maxwellian  case. 
This  is  due  to  a  relative  increase  in  the  0+  density 
for  the  softer  electron  flux  case.  This  may  be  seen 
along  with  overall  behavior  of  the  dominant  ions  in 
Figures  19-22.  The  N0+  density  generally  dominates 
the  total  ion  concentration  but  it  should  be  kept  in 
mind  that  the  density  ratio  of  N0+  to  0^  is  sensitive 
to  the  NO  density  (see  Figure  10). 

We  now  consider  some  of  the  important  UV 
emissions  and  their  properties  for  the  various  inci¬ 
dent  fluxes.  Figure  23  gives  the  calculated  altitude 
profiles  of  the  N2  IN  3914  A  column  emission  rate  for 
viewing  in  the  zenith  direction.  Other  results  to 
follow  also  apply  to  this  observing  direction.  Figure 
24  shows  altitude  profiles  for  01  1356  A.  Substantial 
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<J). (E/p)  el/cm  -s-eV-sr 


300 


IGURE  17.  TOTAL  ION  PRODUCTION  RATES  FOR  THE  INCIDENT  FLUXES  IN  FIGURE  16 


differences  are  seen  between  these  figures.  The  fall-ofl' 
at  low  altitudes  for  01  1356A  is  due  to  photon  absorption 
by  C>2 .  The  overall  decrease  in  its  emission  with  increas¬ 
ing  hardness  in  the  incident  spectrum  comes  from  a  combi¬ 
nation  of  increasing  photon  absorption  by  O2  and  a  rela¬ 
tive  decrease  in  the  available  0  column  density  with 
respect  to  those  of  N2  and  O2 . 

Figures  25-28  show  column  emission  rates  for 
i\'2  3914A ,  01  135GA ,  N2  2P  3371A,  N2  LBH  138-1  A,  and 
h’2  FB11  1325A  for  each  of  the  four  incident  fluxes. 

Here,  we  clearly  see  the  relative  decrease'  in  01  1356A 
emission  relative  to  N2  3914A  as  the  incident  electron 
spectrum  becomes  harder.  The  relative  strengths  of  the 
LB11  bands  also  change,  in  this  case  due  purely  to  dif¬ 
ferences  in  photon  absorption  by  O2.  The  effect  as  seen 
from  a  satellite  above  the  emitting  region  would  be  a  de¬ 
crease  in  the  emission  ratio  of  1381A  to  1325A.  Such 
an  effect  could  be  used  to  help  estimate  the  hardness 
of  the  precipitating  flux.  As  can  be  seen  from  the 
figures,  other  such  emission  ratios  would  also  vary 
which  could  provide  further  information  on  spectral 
hardness . 
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500  eV 


FIGURE  25.  SELECTED  VERTICAL  VIEWING  COLUMN 
RATES  FOR  AN  INCIDENT  MAXWELLIAN 
CHARACTERISTIC  ENERGY  En  =  500  eV 
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1000  ev 
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Appendix  A 


TRANSPORT  EFFECTS  ON  ION  DENSITIES 

We  wish  to  know  under  what  conditions  the  local 
approximation  is  valid  lor  solving  the  time  dependent 
continuity  equation  (Equation  4 ) .  This  is  a  function 
of  the  life  of  the  auroral  event,  species  chemical  life¬ 
time,  altitude,  lateral  extent  of  the  precipitation  pat¬ 
tern  and  strengths  of  an  E  field  and  wind  which  may  be 
present.  We  can  only  speculate  about  the  effect  of  some 
of  these  parameters  since  they  may  be  poorly  specified 
in  a  given  situation.  Nevertheless,  we  can  state  with 
some  confidence  that  the  local  approximation  is  generally 
good  throughout  the  auroral  E-region.  It  will  eventually 
break  down  with  increasing  altitude  due  to  decreasing 
collision  frequencies  (increasing  conductivities)  which 
increase  particle  mobilities.  It  is  to  this  part  of  the 
ion  transport  problem  we  will  direct  our  attention  below. 

The  generalization  of  Equation  (4)  with  trans¬ 
port.  is 

Dr^u.-t)  rVu.i)  -  Ljtl.-O  -V-Jj/e-  (Al) 

dt 

where-  Jj  is  current  density  and  Jj/e  is  flux  in  particles 
cri)2-s  of  either  the  .j^h  type  of  ion  or  ambient  electrons. 
For  the  various  transport,  processes ,  we  wish  to  know  when 
the-  v  •  Jj  term  becomes  comparable  in  magnitude  to 
either  P;  or  L;.  Schunk  and  Walker  (l‘)71)  investigated 


p 


this  problem  and  here  we  shall  follow  their  example  (see 
also  Vallance  Jones  (1974,  Chapter  9)).  The  current  J 
comes  from  the  next  higher  moment  of  the  Boltzmann  equa¬ 
tion,  namely  from  the  momentum  equation.  Assuming  a 
steady  state,  we  have 


-VPJ*rbej 

with  terms 


P  . 
J 


2d 


vj 


f  fc  4  u.  x  §  1  +  n  e.-  v  ■  x  8 

*-  p  .)  j  j  -j  — 

co 

defined  as  follows: 

ion  pressure  of  species 

center  of  mass  velocity  of  the 
ion  drift  velocity  relative  to 
ion  collision  frequency 


3  TVIWP.  V  (A2) 
J  J  J 


atmosphere 

u 


For  electrons,  we  have 


♦  U.  * 


'<V  Vt 


x  B 


Ve 


(A3) 


Fluxes  J^/e  anc*  ^e/e  can  explicitly  represented  in 
these  equations  through  their  definitions: 

( A4 ) 

(A5) 


^  /ej 


ti 


and 


/e. 


0e  Ve 


where  the  charge  terms  e^  and  e  are  positive. 

V.7e  shall  now  restrict  our  attention  to  the 

j  .  £ 

field  aligned  flux  component  - --  where  k  points  up¬ 
ward  along  B  and  to  singly  charged  ions  such  that 

e.  =  e.  For  ions,  we  have 
J 
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A 

*  - 


vA 


^  A  r  \  A 

V?'.  \e  +  wc.  t  -  k  4  fV'fW  \j  •  k 

U  J  «J  sj  — 


] 


Wc  nolo  that  here  four  effects  are  responsible  for  the 
flow  of  ions:  gravity,  diffusion  through  VPj ,  an  E 
field  force  and  mass  motion  by  wind  given  implicitly 
through  u.  The  parallel  component  of  E  may  contain 
an  external  parL  along  with  that  caused  by  the  greater 
field  aligned  diffusion  of  electrons  compared  to  ions. 
This  comes  from  the  electron  momentum  equation  and  has 
t  ho  f  orin 


k  =  -  [  -V?.  fe  -  Ttc.-W 
L  - 


The  external  part  of  the  field  is  responsible  for  cur¬ 
rent  Je.  Substituting  Equation  A7  into  A6,  we  have 


'V 

*  w. 


7,  l  ‘VP  \  -  ^0  V?6-l  -  no^VedfctJ.  4Y)(Mfr.k 
Vo1-  J  n,  7  '  J  4‘ 


We  wish  to  evaluate  these  terms,  their  divergences,  and 
th*-n  compare  the*  total  with  the  chemical  loss  rate  at 
various  altitudes.  To  do  so,  we  must  specify  collision 

k,  u  and  a  model  ionosphere 
and  u  to  be  zero  but  will  note  when 


requenoies  and  v0,  Je 


We  will  assume  J( 
t.tie.v  become  important.  The  following  simple  forms  will 
be  assigned  to  tlu>  frequencies  (see  Jones  (1974) 
(Equation  0.2g)  and  references  therein): 


-H 

o  •-  a  •  to  mo  nn  / 


Vu 


and 


_io 


■pt  -  5  M5 


4  s.^-i  md 


L*v  /  i 

0 
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(A6) 


(A7) 


(A8) 


( A9 ) 


(A10) 


where  nn  is  the  total  neutral  density  (cm  3)  and  pn ;  is 
JijiUn  -J 

reduced  mass  - —± - . 

m.  +  m 
J  n 

For  a  model  ionosphere,  we  shall  use  the  infor¬ 
mation  provided  by  Table  9  and  Figures  9  and  10.  Thus, 
the  electron  and  ion  densities  are  those  calculated  for 
the  incident  electron  flux  shown  in  Figure;  8.  As  noted 
before,  its  energy  dependence  above  *v.  1000  eV  is  similar 
to  that  of  a  2000  eV  Maxwellian  and  ’ts  power  content  is 

O 

^  2  ergs/cm  -s.  The  ion  densities  will  be  taken  from 
the  middle  panel  in  Figure  10. 

The  pressure  terms  will  be  evaluated  in  terms 
of  densities  and  temperatures  as  follows: 

^•Y,  -  W7  t)n)  4  W  fi  ■  (All) 

IPl  J  y-i 

and 

S/'p.W  -  feX  C )jn<  4  WlTL  VJc  ( A12 ) 

6  e  Tz  d* 

J  .  •  k 

Tables  A1  and  A2  give  the  flux  e -  and  con¬ 

tributing  terms  from  Equation  AS  for  N0+  and  0+.  The  ion 
0+  is  included  since  it  is  the  dominant  ion  species  above 
'v  200  km.  As  noted  before,  these  fluxes  do  not  contain 
contributions  from  an  external  electric  field  (Je  =  0)  or 
wind  (u  =  0).  The  calculated  fluxes  are  negative  which 
means  they  are  moving  in  the  downward  direction.  For 
N0+,  this  happens  in  spite  of  the  upward  force  coming 
from  the  negative  slope  of  the  N0+  density.  This  force 
is  overcome  by  the  combination  of  gravity  and  the  ambi- 
polar  force  due  to  downward  movement  of  electrons. 
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TABLE  A1.  N0+  PARTICLE  FLUX  (/cm2-S)  AND 

TERMS  NEEDED  TO  EVALUATE  IT 

O 

- (dynes/cm  ) - 

(s'1)  (cm'2  -  s"1) 


Z(km) 

-VPi  •  k 

— J-VP  •  k 
ne  e 

_  A 

nimig  •  k 

*  R/e 

150 

3.8(1) 

-2 . 3( -16 ) 

-9 . 2 ( - 16 ) 

-2 . 2 ( - 15 ) 

-2.0(6) 

180 

1.1(1) 

2.1(-16) 

-2 . 4  ( - 1 5 ) 

-1.6(-15) 

-9.1(6) 

200 

5.7(0) 

4 . 9 ( - 16 ) 

-1.5(-15) 

-1 . 6 ( - 1 5 ) 

-9.1(6) 

220 

3.4(0) 

8. 1  ( -16) 

-8 . 4 ( - 16 ) 

-1 . 0( -15 ) 

-6.1(6) 

250 

1.6(0) 

4. 3( -16) 

-4. 3 ( - 1 6 ) 

-3. 7 ( -16) 

-4.8(6) 

8G 


Table  AH  provides  a  comparison  between  volu¬ 
metric  rate  changes  due  to  transport  (v  •  Jj/e)  and 
chemical  loss  for  N0+  and  0+ .  We  observe  that  transport 
effects  are  unimportant  to  beyond  250  km  for  the  partic¬ 
ular  aurora  considered.  This  is  generally  consistent 
with  findings  of  Schunk  and  Walker  (1971)  although  they 
considered  a  somewhat  softer  incident  electron  spectrum 
containing  more  energy. 

To  complete  this  discussion,  we  wish  to  indi¬ 
cate  at  what  strengths  flow  due  to  an  external  electric 
field  and  wind  become  comparable  with  flow  due  to  terms 

already  considered.  From  Equation  A8 ,  the  field  term  is 
n  j  J  e  •  k 

given  by  ^  m(>  v(.  - .  It  must  have  values  on  the 

order  of  2  x  10-^^  and  5  x  10-1^  dynes/cm^  between  200 
and  250  km  respectively  for  N0+  and  0+  to  be  comparable 
to  the  sum  of  the  other  terms  (see  Tables  A1  and  A2). 

J  «  ^ 

These  values  correspond  to  electron  currents  (  ■£-~ - )  in 

o  e 

the  range  of  100  pA/m  . 

To  estimate  at  what  strength  the  parallel  com¬ 
ponent  of  the  wind  becomes  important  for  the  given  set 
of  auroral  conditions  under  discussion,  we  simply  need 
to  examine  flow  velocities  corresponding  to  the  currents 
in  Tables  A1  and  A2.  They  are  roughly  between  5  and 
20  m/s  in  the  altitude  range  from  200  to  250  km.  Thus, 
vertical  wind  speeds  on  the  order  of  10  m/s  are  called 
for  to  impact  on  the  analysis  presented  in  this  appendix. 
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TABLE  A3.  COMPARISONS  OF  TRANSPORT  AND  CHEMICAL 
LOSS  RATES  FOR  NO+  AND  0+ 


Z(km) 

V  •  JN0+/e 

anenN0+ 

V  •  JQ+/e 

( k  jno2+^2nN2 )  r'o+ 

150 

-1.5(0) 

9.7(2) 

-- 

2.2(2) 

180 

-1.7(0) 

3.2(2) 

-3.5(0) 

9.1(1) 

200 

0.0 

1.7(2) 

-8.8(0) 

5.5(1) 

220 

7. 5( -1 ) 

1.1(2) 

-2.2(1) 

3.1(1) 

250 

1.0(0) 

4.2(1) 

-7.0(1) 

9.9(0) 
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